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Abstract 

Atomic ensembles have many potential applications in quantum information science. Owing to 
collective enhancement, working with ensembles at high densities increases the overall efficiency 
of quantum operations, but at the same time also increases the collision rate and leads to faster 
decoherence. Here we report on experiments with optically trapped Rh atoms demonstrating 
a 20-fold increase of the coherence time when a dynamical decoupling sequence with more than 
200 pi-pulses is applied. We perform quantum process tomography and demonstrate that using 
the decoupling scheme a dense ensemble with an optical depth of 230 can be used as an atomic 
memory with coherence times exceeding 3 sec. 



Cold atomic ensembles can be used as an interface between matter and photonic qubits 
in quantum networks [1, 2], and in recent years vast experimental advances in this direction 
have been reported [3-9]. Quantum information, which is mapped into the coherence be- 
tween two atomic internal states, is gradually lost due to inhomogcneities and fluctuations 
in the energy difference between these states. For trapped atoms the inhomogcneities are 
caused by differential light shift in optical traps [10] or by differential Zeeman shifts in mag- 
netic traps, and by mean- field density dependent interaction shifts [11]. Fluctuations arise 
due to collisions which are inherent to the high densities required to achieve a good overall 
efficiency of quantum operations [12, 13]. Though fiuctuations at low frequencies can be 
overcome by a single population inverting pulse - the celebrated coherence echo technique 
[14, 15] , as the collision rate increases this is no longer possible due to higher frequency com- 
ponents. Dynamical decoupling theories generalize this technique to multi-pulse sequences 
by harnessing symmetry properties of the coupling Hamiltonian [16-20]. Though dynamical 
decoupling was demonstrated in several experiments [21-26], its use with atomic ensembles 
remains unexplored. In addition to its practical importance, this exploration is of theoretical 
interest since in trapped atomic ensembles the energy distribution is non-Gaussian and the 
fluctuations originate in self-interactions and not in a noisy external environment. 

We consider atoms trapped in a conservative optical potential. The effective single par- 
ticle Hamiltonian for atoms with internal states designated by |1) and |2) is given by 

H = h[uJo + S{t)]\2){2\+hn{t)\2){l\ + h.c. , (1) 

where coo is the free space transition frequency between the states, 6{t) is a random fre- 
quency detuning sequence whose nature is determined by the potential inhomogcneities and 
collisions, and is the external control field which is used for the dynamical decoupling. 
Starting with an initial state |'0(O)) = 2~^/^(|l)-|-|2)) and no external control fields, the wave- 
function at any given time is given in the rotating frame by \ip{t)) = 2~^/^(|l) + e"""^*^*^ |2)), 
where the phase difference is given by (pit) = S(t)dt. A schematic plot of three real- 
izations of 0(t) is given in Fig. la (top), and it can be seen that the phase difference is 
accumulated in a constant rate between collisions [27]. The ensemble coherence is charac- 
terized by the function C{t) — where pi2 is the off-diagonal element of the reduced 
two-level density matrix [20]. As an example, for a Gaussian phase distribution, P^, with 
a standard deviation cr,/, we obtain C{t) = e , which shows that the coherence decays 
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as the width of the phase distribution increases. The effect of a population inverting pulse 
(tt- pulse) is to change the sign of 5, and a train of such pulses lead to a much narrower phase 
distribution and slower decoherence, as depicted in Fig. la (bottom). 

The experiments are carried out with cold ^^Rb atoms trapped in a far-off-resonance 
laser (see Fig. lb, and more details in the supplementary information). The two relevant 
internal states are |1) — \F — l;m/ — 1) and |2) — \F — 2;m/ = 1) in the 5^5'i/2 manifold, 
which are, to first order, Zeeman insensitive to magnetic fluctuations in the applied magnetic 
field of 3.2G [11]. The external control is done by means of two-photons transition (MW- 
RF photons), and the detection is state sensitive [28]. By gradually lowering the trapping 
laser intensity we reach the experimental conditions, at which we have 275, 000 atoms at a 
temperature of l.T/iK, phase space density of 0.04 and an average collision rate of 100s~^. 
The typical inhomogeneous decay time as measured in a Ramsey-like experiment is ~ 150ms. 
The peak optical depth for a non-polarized resonant light is ~ 230. 

We employ a Carr-Purcell-Meiboom-Gill (CPMG) decoupling scheme [29] and show in 
what follows that for collisional detuning fluctuations it is virtually optimal. The pulse 
sequence is composed of n 7r-pulses at times t/. = ^|^^ where k = 1 . . .n (see Fig. Ic), 
and we characterize it by the effective frequency fon — We study the effect of the 
dynamical decoupling scheme by performing a quantum process tomography (QPT) [30]. 
QPT enables us to reconstruct the x-matrix which gives a convenient way to calculate 
the density matrix after the process, Pout, in terms of the initial density matrix, pin, by 
Pout = £[pin] = Y.k,i^kpinEjxkh whcrc E = {i,X,-iY,Z) with {i,X,Y,Z) being the 
Pauli matrices. In the experiment we start with a set of initial states after which we apply 
the decoupling scheme and measure pout by quantum state tomography (for more details see 
the supplementary information) . The results of a QPT of a dynamical decoupling sequence 
with foD — 35Hz is depicted in Fig. 2a. There are two distinctive decay timescales for the 
equatorial plane and the 2;-axis, which corresponds to phase damping noise processes and 
depolarizing noise processes (Ti), respectively. The former originates from the fluctuations 
in 5 and it is the dominant noise process which determines the ensemble coherence time, Tc, 
and is quantified by C{t). Depolarization process is induced by inelastic collisions [31], and 
its typical timescale is measured to be Ti = 6s. The worst case fidelity of the ensemble as a 
quantum memory, defined as = min('0|£^ [j'^) (■0|] j'^), is calculated from the measured x~ 
matrix to be = 0.83, 0.74 and 0.64 for 1, 2 and 3 seconds, respectively, which corresponds 
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FIG. 1. The experimental setup and dynamical decoupling sequence, a, Schematic drawing of the 
evolution of the relative phase between the two atomic internal states with and without the dynam- 
ical decoupling pulse sequence. Without dynamical decoupling (top), after a collision the average 
potential energy of the atom is changed and therefore also its average rate of phase accumulation, 
6. The inset shows the small oscillations due to the fast periodic atomic motion in the trap. Since 
the oscillation period is shorter than all relevant timescales in our experiment, we shall consider 
only 5{t) which is averaged over several oscillations [27]. The dynamical decoupling pulse sequence 
is plotted in the lower graph, where each pulse is a vr pulse. With dynamical decoupling pulse se- 
quence (plotted in the bottom), the spread of the phases is much smaller (note the different graph 
scales), b, The experimental apparatus. We laser cool ^'^Rb atoms and trap them in a crossed 
red-detuned laser beams configuration. We employ state sensitive detection using a detection beam 
and a photomultiplier tube (PMT), and measure the density and temperature using absorption 
imaging on a CCD camera, c. The experimental pulse sequence starts with a state preparation, 
followed by a train of population inverting pulses with alternated phases {vr, vr, — vr, — vr, ...} to min- 
imize the accumulation of errors due to pulse width and frequency inaccuracies, and a final state 
detection. The duration of each vr pulse is ~ 0.5ms and its average fidelity is ~ 0.995. 



to an exponential decay timescale of = 2Asec. The contraction of the Bloch sphere is 
symmetric in the equatorial plane, which indicates that the decoupling scheme is insensitive 
to the stored superposition. We demonstrate this point with a direct measurement in which 
we start with two orthogonal initial states in the equatorial plain: \ipi) = :^(|1) + |2)) and 
\ip2) = + ^^'^^'^ For both states we scan the phase of a final tc/2 pulse added to 
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the sequence and measure the population at |2) normahzed to its initial value. The results 
depicted in Fig. 2b exhibit the same contrast and preserve the 7r/2 phase shift between the 
two initial states. 

The decay of the coherence with the dynamical decoupling pulse sequence and assuming 
a Gaussian phase distribution can be obtained in a system-reservoir framework [20, 32] : 

C(t) = e"-^o°°'^^«H^M)A'^^ , (2) 

where the argument is the overlap integral between the fluctuations power spectrum, Ss{uj) = 
J^(5(t)5(0))e*'^*(it, and a filter function which encapsulates the information on the dynam- 
ical decoupling pulse sequence and is given by F[ujt) = ^| ^^=o(~-'-)^(^*'^*''^^ ~ g*'<jifc^|2 ^j^^j^ 
to — and tn+i — t. For each atom 5{t) is a sequence of constant detunings connected by 
"jumps" which occur after each coUision. Since the collision times follow Poisson statistics, 
the detuning correlation function decay exponentially ^5{t) — {6{t + 0)6(0)) — (7|e~'"l*l, 
where F"^ is the correlation time of the detuning and as is the standard deviation of the 
detunings distribution. The power spectrum is given by Ss{io) = J^^^s{t)e'^^^dt = ^2^^2 ■ 
By solving numerically Eq. (2) and leaving the {ti}^^i as free parameters, we find that 
the optimal decoupling sequence for a Lorentzian power spectrum is given by U — ^^l^l^ t, 
where i — l...n and 0.5 < 77 < 1 is a numerical factor which depends on n and t. For ^ <^ 1 
we find r] ~ 0.5, for which we retrieve the CPMG pulse sequence. Furthermore, even when 
^ ^ 1 the coherence time with the CPMG pulse sequence differs by less than 1% from the 
optimal value. We have tested theoretically and experimentally other dynamical decoupling 
schemes, and in particular the one suggested in Ref. [19], and verified that they are indeed 
inferior to the CPMG sequence in our Lorentzian fiuctuations power spectrum (for more 
details see the supplementary information). 

We measure C{t) directly by the following sequence: a short 7r/2 pulse with a phase </? 
prepares the atoms in the superposition |^) = + e*'^ |2)), following are the decoupling 

pulses and finally the coherence is measured from the length of the Bloch vector measured 
with a quantum state tomography. Though C{t) docs not have to follow, a priori, some well 
defined function, the experimental results depicted in Fig. 3a show that the data is well 
fitted by an exponentially decaying function e~*/'^'=. The exponential decay is expected in the 
Markovian limit, where the decay timescale is much larger than the fluctuations correlation 
time [20]. This is indeed the case in our experiment, as can be seem in the inset of Fig. 3a 
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FIG. 2. Quantum process tomography (QPT) of dynamical decoupling in cold atomic ensemble, a, 
We perform QPT of dynamical decoupling with /d_d = 35Hz and reconstruct the x-matrix which 
defines the process £[p\. Any single-qubit density matrix p can be mapped to a point in space 
r{p) = {tr{X p),tr{Y p),tr(Zp)), with {X,Y,Z) being the Pauli matrices. The colors and lines are 
chosen for the initial states, pm, which lie on a sphere with a radius 1. For each of these states 
we calculate the process outcome, pout = £[pin], and plot it with its initial color at f{pout)- The 
contraction of the Bloch sphere is more pronounced on the equatorial plane which shows that the 
main noise process is phase damping. Slower depolarizing process (Ti) which uniformly shrinks the 
sphere are also apparent. This can be also seen in higher values of Xzz compared to Xxx and Xyy 
The fidelity of the ensemble as a quantum memory is calculated to be 0.83, 0.74 and 0.64 after 1, 
2 and 3 seconds, respectively, yielding a coherence time of Tc = 2.4 sec. There is also a rotation of 
the sphere around the |1)-|2) axis at a rate of ~ 9°sec~^ due to small inaccuracies in the control 
field, b. Storage of two orthogonal initial states in the equatorial plane: l^/^i) = + and 

14^2) = + We apply the decoupling scheme and add to it another pulse independent 

of the initial state with a phase and duration that are chosen to correct for the small rotation of the 
Bloch sphere. We measure the population at |2) after 3 sec, normalized to the initial population, 
versus the phase of a tt/2 detection pulse. The phase difference between the two initial states is 
conserved and also the contrast which corresponds to a higher Tc = 3 sec due to the correction of 
the rotation. The fringe contrast is not centered to 0.5 due to inelastic m-changing transitions in 
the hyperfine level F = 2. 
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FIG. 3. The ensemble coherence time with a CPMG dynamical decoupling pulse sequence, a, 
The coherence left in the ensemble without dynamical decoupling {foD = OHz) and with two 
representing dynamical decoupling pulse rates after different waiting times, normalized to the 
initial stored coherence. From the fit to the exponential decay function (dotted line) we extract 
the coherence time. The inset shows calculations according to Eq.(2) of the coherence time in 
a semilog scale, and demonstrates the expected exponential decay behavior, b. The parabolic 
dependence of the coherence time on the dynamical decoupling pulse rate. The experimental data 
(blue circles) agrees well both with the theoretical prediction of Eq.(2) taking a Lorentzian power 
spectrum (solid line) with as = 23.8sec~^ and F = 37.5sec~^ which were extracted from measured 
quantities and with no fit parameters, and with molecular dynamics Monte-Carlo simulations 
done with 1000 atoms in a 3D harmonic trap with an inhomogeneous decay time similar to the 
experiment (red squares). The errorbar is a estimated from the fits shown in (a), c, The dephasing 
rate is measured for different collision rates, demonstrating the linear dependence of Eq.(3). These 
experiments were done with /dd = 8Hz. 

where we calculate C{t) using Eq. (2) with several values of Jdd] The lines on a semilog 
scale show deviations from linear fits only for t ^ F^^. We also measure the dependance of 
the coherence time on the dynamical decoupling pulse rate and the results are shown in Fig. 
3b. We observe a quadratic increase of the coherence time versus /dd up to 35Hz, where 
we obtain a 20-fold improvement with coherence times exceeding 3 sec. 

In order to explain these results we present a qualitative model for the coherence time. 
Without collisions the inhomogeneous dephasing rate is proportional to as. For simplicity we 
assume that if a collision did not occur between two consecutive vr-pulses the inhomogeneous 
broadening is averaged out. If a collision occurred, however, the width of the ensemble 
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phase distribution increases by ~ foo'^s- The number of coUisions up to a time t is Vcoit^ 
and since we add random variables (i.e. the accumulated phase), the width of the phase 
distribution increases as a square root of time: A$(t) ~ /^^cr^vT^- For cold coUisions in 

3D harmonic trap, the relation between the collision rate and the relaxation rate was shown 
to be [33] Tcoi = 2.7 - T. The coherence time, Tc, is the time for which the width of phase 
distribution is on the order of 1, yielding 

Tc-foD'yf'^-' , (3) 

with a parabolic dependence on Jdd- This result can be also obtained from Eq.(2) by 
approximating ~ Soiraci'^t — 2nfDDt) and using the Lorentzian power spectrum. 

Exact calculations of Tc using Eq.(2) without fitting parameters are presented in Fig. 3b 
in good agreement with the experimental data. The calculations are done with a Lorentzian 
power spectrum where the parameter F and as are calculated from measured quantities. 
F is extracted from the collision rate which is calculated using the measured temperature, 
number of atoms and trap oscillation frequencies. The parameter is measured in a Ram- 
sey experiment at very low densities, where the collisions can be disregarded and as can 
be extracted from the measured dephasing rate (for more details see the supplementary 
information). We also perform Monte-Carlo simulations, where we solve for the classical 
motion of atoms in the trap including collisions, and calculate the Ramsey signal by tracing 
each atom's accumulated phase along its trajectory. The results of the simulations are also 
depicted in Fig. 3b, and agree well with both theory and experiments. We conclude that the 
effect of collisions can be indeed formulated as an effective single spin Hamiltonian coupled 
to a reservoir. Moreover, although the detunings of atoms trapped in a 3D harmonic trap 
are not normally distributed (for more details see the supplementary material), the distri- 
bution of their accumulated phase can be well approximated by a Gaussian owing to the 
validity of the central limit theorem and large number of collisions involved. 

Another outcome of Eq.(3) is the linear dependence of the coherence time on F^^. In the 
experiment we change F by reducing the density and collision rate while keeping the tem- 
perature, and therefore as, unchanged [27]. This is accomplished by reducing the intensity 
of the cooling lasers in the trap loading phase. In Fig. 3c we plot versus the average 
collision rate for a pulse rate oi foD — 8Hz. As expected, the coherence time is inversely 
proportional to the collision rate. 
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In conclusion, we have demonstrated that dynamical decoupling can substantially in- 
crease the coherence time of a dense optically trapped atomic ensemble. In the current 
work the ensemble was treated as an effective single spin system which accounts for storage 
schemes based on collinear EIT [2]. A natural extension is the application of dynamical 
decoupling to Raman scattering schemes in which a global coherence is created between the 
atoms [1, 3]. Another promising prospect lies in novel hybrid approaches to quantum compu- 
tation combining atomic ensembles and superconducting devices [34], where the application 
of dynamical decoupling could reduce the error probability during the characteristic 100/xs 
single qubit gate to less than 10~^ - below the current estimated threshold for a fault-tolerant 
quantum computation. 

We thank Roee Ozeri, Guy Bensky, Goren Gordon, Gershon Kurizki and Rami Pugatch 
for helpful discussions. We acknowledge the financial support of MIDAS, MINERVA, ISF, 
and DIP. 
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THE EXPERIMENTAL SETUP 



Initially ~ 10^ atoms are trapped and cooled in a magneto-optical trap, and further 
cooled by Sisyphus and Raman sideband techniques. We use rapid adiabatic passage 
with constant RF radiation and a ramped magnetic field to transfer the atoms from 
state |5^S'i/2,-F = l;m/ = l) ending with 80%of the atoms at and the rest in state 
15^5*1/2, -F = 0;m/ = O). The atoms arc loaded into an optical dipole trap created by two 
horizontal crossing beams at an angle of 28°, creating an oval trap with an aspect ratio of 
1:3.9 with a measured radial oscillation frequency of 27r-330Hz. The 50//m waist laser beams 
originate from a single frequency Ytterbium fiber laser at 1064nm. Their polarization is 
parallel to the magnetic field, and frequency differ by 120MHz to eliminate standing waves. 
The external control is composed of RF radiation at 2.15MHz and microwave radiation at 
~ 6. 8GHz, both locked to an atomic standard. Another effect which should be taken into 
account is the levels shift induced by the MW field [1]. We have carefully measured this 
shift, which in the maximum MW power reaches to ~ 50Hz, and we take it into account 
when setting the frequency of the external fields. We carry out evaporative cooling for 2.5sec 
to a laser power of 0.16W and back to a final value of 0.64W. The spontaneous scattering 
rate is less than ls~^ and the trap lifetime is longer than 5s. We measure the temperature 
by an absorption imaging after a time of ffight. 

QUANTUM STATE TOMOGRAPHY 

In quantum state tomography the super-operator matrix x is determined by implementing 
a complete set of qubit input states. For a single qubit this is conveniently done by a 
measurement of tr[Xp\, tr[Yp] and tr[Zp\ [2]. We measure the population at state |2) with 
and without another tt pulse, and the difference of the two values normalized to the initial 
population at state |1) gives MM±1 obtain the other two projections by apply a 7r/2 
pulse in the suitable axis and then measure tr[Zp]. 

DEPHASING IN A 3D HARMONIC TRAP 

The dephasing in a 3D harmonic trap without collisions is given by C (t) 
[3], and solving C(ri) = e"^ we obtain ti = \/3 (e^/^ — l)aj^ fa 1.69 • a^^ 
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FIG. 1. Coherence time versus r] for f£,£, = 35Hz and a lorentzian power spectrum with the typical 
experimental parameters F = 37. 5s^^ and = 23.8,s~^. The calculation of is done at t = Is. 

OPTIMAL DECOUPLING PULSE SEQUENCE FOR A LORENTZIAN SPEC- 
TRAL FUNCTION 

As explained in the text, we want to find the pulse sequence timing series {ti}^^^ which 
maximizes the coherence as calculated by Eq.(2). We use a nonlinear minimization method 
(Nelder-Mead) implemented in MATLAB software for the minimum decoherence. We have 
consistently obtained a time series close to the form given in the text: tj = ^^l^^^ t^ where 
i = l...n and 0.5 < r/ < 1 is a numerical factor which depends on the parameters of the 
problem. For F, the difference of the time series for the optimal rj to the one with 

1] = 0.5 (which corresponds to the CPMG scheme) is less then 1%. As an example we plot 
the coherence time as a function of rj for foD = 35Hz in Fig. 1. The difference in this 
case between Tc of t] = 0.5 and the optimal value is ~ 0.14%. In Fig. 2 we plot the same 
calculation for /dd = lOHz. Here the difference between the optimal value and the value at 
T] = 0.5 is ~ 0.93%. Therefore the value i] = 0.5 was chosen to be used in the experiment. 

In regard to the Uhrig dynamical decoupling scheme (UDD) [4], we plot in Fig. 3 a 
comparison between CPMG and UDD as a function of the number of vr-pulses. As can 
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FIG. 2. Coherence time versus rj for fDD = lOHz and a lorentzian power spectrum with the typical 
experimental parameters F = 37.5s~^ and as = 23.8s~^. The calculation of Tc is done at t = Is. 
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FIG. 3. Coherence time versus the number of 7r=pulses for a CPMG dynamical decoupling scheme 
and the UDD scheme, for a lorentzian power spectrum with the typical experimental parameters 
F = 37.5s~^ and as = 23. 8s^^. The calculation of Tc is done at t = Is. 
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be clearly seen from the graph, the UDD scheme is inferior to the simpler CPMG scheme 
for a Lorentzian bath spectrum. This is not surprising since the UDD is expected to be 
superior in spectral functions that exhibit a cutoff in frequency. This point was already 
discussed in Ref. [5]. We have also tested this experimentally for several number of pulses 
and consistently obtained coherence times with the UDD scheme which are 10%-20% lower 
than the coherence times obtained with the CPMG scheme. 
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